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SYNTHESIS  OF  THE  FORMS  OF  ATMOSPHERICS 
D.  S.  Fllgel* 

Introduction 

Atmospherics  formed  from  lightning  discharges  are  a  convenient 
object  for  Investigating  long  and  super-long  electromagnetic  waves. 

In  previous  works  at  the  Institute  of  Terrestrial  Magnetism^  the 
Ionosphere 4  and  Radiowave  Propagation  [IZMIRAN]  [1,  2]  It  was  shown 
that  the  phase  velocity  of  the  electromagnetic  waves  conqjrlslng  atmos¬ 
pherics,  and  the  effective  conductivity  of  the  lower  Ionosphere,  can 
be  determined  by  analyzing  the  phase  (and  to  some  extent  the  ang)ll- 
tude)  characteristics  of  single  atmospherics.  A  reasonably  correct 
solution  of  the  opposite  problem  was  also  of  great  Interest,  viz., 
to  calculate  theoretically,  on  the  basis  of  a  specific  model  of  the 
Ionosphere  evolving  from  experimental  data  (see,  e.g.,  [5])  and  on 
the  basis  of  experiments  on  the  spectral  properties  of  the  source 
(see  [2,  4,  5])#  the  forms  of  atmospherics  E(t,r)  expected  at  various 
distances.  ConQ}arlson  of  the  signals  thus  calculated  with  the  atmos¬ 
pherics  observed  experimentally  would  make  It  possible  to  verify 
first,  the  correctness  of  the  general  theoretical  assumptions  used 
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for  this,  and  second,  to  explain  whether  the  effective  parameters 
obtained  theoretically  for  the  selected  Ionosphere  model  were  observed 
during  the  experiment,  and  If  so,  how  frequently. 

The  form  of  the  atmospherics  E(t,r)  has  already  been  calculated 
by  Sudden  [6]  for  a  distance  r  =  1000  km  from  the  radiator.  As  the 
radiation  source  he  used  a  dipole  for  which  the  relation  between  the 
electric  moment  and  time  has  the  form  of  a  step,  l.e.,  he  selected  a 
source  which  did  not  reflect  the  form  E(t)  of  the  close  discharges 
5].  It  was  also  assvimed  that  the  Ionosphere  has  Identical  con¬ 
ductivity  over  all  the  frequencies  con^jrlslng  the  atmospherics.  At 
the  same  time,  Al'pert  [3]  has  shown  that  In  the  range  of  frequencies 
comprising  the  atmospherics  spectrum  the  effective  conductivity  of 
the  Ionosphere  causing  the  propagation  of  these  waves  changes  by  more 
than  100- fold.  In  addition.  In  Budden’s  work  [6]  all  calculations 
were  made  on  the  assumption  that  when  r  =  1000  km  only  one  null  mode 
participates.  However,  at  frequencies  f  >  15-20  kc  at  this  distance 
It  Is  necessary  to  consider  no  less  than  2-3  modes  [3].  Thus  although 
the  signal  obtained  by  Budden  [6]  was  generally  similar  in  form  to 
the  atmospherics  observed  experimentally,  we  could  draw  no  further 
conclusions  from  this.* 

Walt  [7]  has  recently  made  similar  calculations.  He  derives  the 
forms  of  the  atmospherics  calculated  for  various  distances  from  a 
source  for  the  correct  form  of  the  source.  Nonetheless,  Walt  [7], 
like  Budden  [6],  assianes  Identical  Ionospheric  conductivity  for  all 
frequencies,  and  the  calculations  are  made  only  for  one  mode.  There¬ 
fore  the  obtained  signals  do  not  show  any  dependence  on  distance  and 
are  \xnsuitable  for  any  rigorous  conparlson  with  experimental  data. 

*  Budden  [6]  did  not  set  this  as  one  of  his  tasks. 


FrD-TT-62-1195-1+2+4 


-2- 


We  learned  of  Walt's  work  [7]  only  after  we  had  finished  and  formu¬ 
lated  our  Investigations;  these  were  reported  on  briefly  [8,  9]  before 
Walt's  work  was  published.  The  present  article  contains  more  detailed 
results  of  our  work. 

We  calculated  the  atmospherics  for  various  distance  from  a  source, 
taking  Into  account  the  necessary  number  of  modes  and  the  dependence 
of  effective  conductivity  of  the  ionosphere  on  frequency.  We  give 
an  analysis  of  the  calculated  and  observed  forms  of  the  signals  from 
a  number  of  their  external  characteristics,  and  a  more  profound  analy¬ 
sis  Is  given  by  con^arlson  of  their  phase  characteristics.  We  have 
confirmed  the  possibility  [l]  of  determining  In  each  specific  case, 
by  such  a  con^arlson,  the  effective  Ionospheric  conductivity  causing 
the  propagation  of  the  respective  waves.  The  realization  of  this 
method  can  apparently  be  vised  for  practical  purposes.  However,  Its 
accviracy  can  be  evaluated  only  on  the  basis  of  direct  measvirements  of 
the  same  effective  Ionospheric  parameters, 

1.  Calculation  Method 


To  calculate  the  foiro  of  the  signal  received  at  each  disteince 
r,  l.e,,  to  determine  E  vs.  t  and  r,  we  must  calculate  the  Fourier 
Integral 


E{t,  r) 


(1) 


where 


g  («,.  r)  =  B  («.  r)  e-H--  ^  (krSn)  P  (C„)  (  2  ) 

n  9 

Is  the  electric  or  magnetic  field  of  a  wave  of  frequency  <0  above  the 
earth's  surface  [5],  while 

Qo  (<".  0)  ==  Ag  (u),  0)  •>  (  5  ) 


-> 
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Is  the  spectral  radiation  density  of  the  source  in  its  immediate 
vicinity. 

It  has  been  shown  [4,  5]  that  the  form  of  the  signals  excited 
by  lightning  in  the  immediate  vicinity  of  the  discharge  is  described 
quite  well  by  an  expression  of  the  form 


£o(/.0)  =  e— '  -e-P', 

while  from  literature  data  on  the  times  of  the  rise  and  fall  of  a 
signal^  Borodina  [2]  obtained  the  values 

a  =  (0.5— p  =  (0.5-20).10« 

and  showed  that  a  change  in  a  and  ^  within  the  indicated  limits  has 
very  little  effect  on  the  anqjlltude  characteristic  of  the  spectrum, 
although  its  phase  characteristic  changes  considerably.  Nonetheless, 
we  can  select  a  "standard"  source  characterized  by  the  values  a  =  10^ 
and  p  *  10-^  whose  phase  characteristics  are 


f.(u).O) 


_ _ 

-  —  (P  +  «> 

~  ajj  —  U.»  ■ 


(4) 


As  for  g(a),r),  the  Kankel  function  in  e^qjression  (2)  can  be 
replaced  by  its  asyTi?)totlc  representation,  when  |krS^|  >  4-5.  If> 
in  addition,  the  field  strength  is  expressed  in  mv/m,  the  current 
I  in  amperes,  and  h  and  r  in  kilometers,  expression  (2)  assumes  the 
form 


g  ((i>,  r)  =  B  (u),  r)  - 


(5) 


here  the 
where 


"complex"  sine  of  the  angle  of  incidence  s„  -- 
is  the  root  of  the  equation  of  the  poles 

2hkC„  !  i  In  R  =  2n'r, 

0C„  IHkR  (C„) 


=  -a 


(5a) 

(5b) 
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R  is  the  Fresnel  reflection  factor,  a  function  of  the  complex  Ionos¬ 
pheric  refraction  factor  and  h  Is  the  height  of  the  Ionosphere; 
k  *  (o/c;  IZq  Is  the  product  of  the  source  current.  In  amps,  and  the 
operating  height  of  the  radiator;  and  X  Is  the  wavelength  In  kilo¬ 
meters  . 

For  all  practical  purposes,  when  f  >  500  km,  expression  (5)  Is 
suitable  for  calculating  g(a>,r).  In  all  other  cases,  precise  values 
of  the  Hankel  functions  must  be  used  [10]. 

Since  the  aspect  of  the  function  E(t,r)  calculated  using  (l)  Is 
not  a  function  of  the  linear  terms  of  the  arguments  of  the  Integrands, 
to  calculate  E(t,r)  we  must  use  the  Incomplete  phase  <l'((D,r),  which 
has  the  form 

4*  (m,  r)  =  krSn  -  arg  P  (C„) - ^  . 

and  the  auxiliary  phase,  defined  by  the  expression 

Therefore,  In  the  calculations  we  will  deal  only  with  f{(a,r), 
Thvis,  Integral  (l)  can  be  represented  In  the  form 

+00 

+00 

+  <  gir  J  ®  (">  0)  M  +  ?) 


where 

<P  =  'Po(<*)>0)-V'(c»>,r). 

Now  we  can  show  that  the  product  B(aj,r)AQ(a>,0)  Is  an  even  func¬ 
tion  of  the  frequency,  so  that 

+«■ 

/  B(aj,r)A  {a),0)Bln(a)t  +  9)da)  =  0 
*00  ^ 
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and  the  problem  reduces  to  one  of  calculating  the  Integral 
E(t,r)  =  —  /  B(a),r)A  (a),0)cos(ajt-q))dtD. 

"A 


To  calculate  Integral  (7)  It  Is  necessary.  In  principle,  to  know 
the  values  of  the  Integrand  from  zero  to  Infinity,  However,  we 
selected  Integration  limits  from  f  =  5  to  f  =  10-^  cps.  The  upper 
limit  corresponds  to  the  frequency  beginning  with  which  the  module 
of  the  Integrand  assvimes  small  values  and  rapidly  decreases.  The 
lower  limit,  generally  speaking,  can  be  reduced  to  zero.  However, 
the  experimental  data  accviraulated  during  Investigations  of  close 
thunderstorms  were  obtained  using  a  device  which  had  a  pass-bsind 
beginning  at  100  cps;  therefore  we  do  not  know  how  accurately  (4) 
approximates  the  lightning  discharge  at  f  ^  100  cps.  If  we  also  take 
Into  consideration  the  fact  that  the  lowest  frequencies  for  which  we 
Investigated  atmospherics  are  of  the  order  of  several  cps.  It  becomes 
clear  that  It  Is  quite  sufficient  to  limit  ourselves  to  the  selected 
value  of  the  lower  limit  f  *  5  cps. 

The  wave  numbers  kS^j  for  the  Integration  Interval  5-10^  cps  were 
calculated  as  follows.  For  f  =  500-50,000  cps  we  used  the  values  of 
obtained  by  Al’pert  [5I  for  a  number  of  discrete  frequencies,  while 
for  f  <  500  cps  emd  f  >  50,000  cps  the  values  of  were  calculated 
using  approximate  formulas  derived  from  (5a),  When  5  cps  <  f  <  500 
cps  the  following  formulas  give  a  good  approximation: 


cl^ 


2  y  z  hk  I 

Cn  =  «*  -1  2Cj,  .S'ni 

S„.=5r-  V?^T(I 


hJi  I  ,5 


1 


\  1) /I*  Kfi 

_ i_ 

^  2  (i--\)iik 


(8) 

(9) 
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where 


n  =  1 ,  2...,  Sn  ■Sni  +  iSntt 
ft  — 


a 


nn 
hk  ' 


h  Is  the  height  of  the  Ionosphere,  In  km,  and  k  =  co/c.  Formulas  (8) 
were  obtained  under  the  asstunptlon  that  the  Ionospheric  conductivity 
a  has  a  rather  high,  but  finite,  value,  l.e.,  that  Wa/eo)  »  1,  where 
e  Is  the  relative  permittivity.  When  f  ^  5  cps 


K  2  V'^hk  KP  ’ 


(10) 


The  reasonable  accviracy  of  Formulas  (8)  can  be  seen,  e.g.,  from 
a  coii?)arlson,  for  f  =  500  cps,  of  the  values  of  and  obtained 
using  (8),  and  from  the  acciirate  formulas  (see  [5]).  Although  the 
criterion  for  the  feasibility  of  the  above-derived  formulas  Is  the 
Inequality 


and  In  the  case  given  below  (see  Table  l)  this  yields  3.75  >  there 
Is  nevertheless  good  coincidence  among  the  corresponding  data. 


TABLE  1 


c, 

c, 

5.. 

-S.. 

Calculation  trm  aaourato  feiwlas 
(3) 

0,15 

0,343 

1,0476 

0,0491 

Calculation  fMm  appracloato 
fomilas  (8) 

^  o.m 

0,341 

1 ,0481 

0,0481 

As  for  the  calculation  of  for  f  >  30  In  this  case  to 
solve  (5a)  we  can  use  (when  satisfying  the  specific  criteria  [5])  the 
expression  for  the  coefficient  of  reflection  of  a  plane  wave  from  a 
nonuniform  layer  recommended  by  Al'pert  [l]  for  C  «  1,  viz.. 
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(11) 


For  this  we  calculated  the  dependence  of  the  coefficient  of 
reflection  on  C  and  the  frequency  for  the  selected  Ionosphere  model 
using  the  Epshteyn  formula  for  the  transition  zone  [?].  Then  the 
complex  values  of  R(C)  were  approximated  using  expression  (11 )j  this 
made  It  possible  to  determine  and  Bg.  Then,  from  the  equation  of 
the  poles  (5a)  we  found  and  Cg,  for  which  we  obtained  the  equations 


2hkCin  —  2Hit  —  n  —  (BjCjn  4-  BtCyt),  ( 12  ) 

B,C,„-(Bj+  2/iA>)C„=0. 

As  a  result  we  obtained  the  desired  values  of  for  the  entire 
range  of  frequencies  f  =  5-10^  cps.  Table  2  gives  the  values  of 
and  P(Cq)  for  the  null  modej  the  corresponding  values  of  and  P(Cj^) 
for  f  »  500-30,000  cps,  n  =  1,  2,  etc.,  have  been  guven  by  Al'pert 
[3].  Plgvires  1-8  show  graphs  of  S^g,  |P(C)|,  and  arg|p(C)|  vs. 
the  relationship  between  the  electron  concentration  and  the  number 
of  collisions  N/v,  obtained  by  Al»pert  [3l  for  the  null  mode  by  solv¬ 
ing  the  pole  equations  (5a)  for  the  Index  of  Ionospheric  refraction 


n* 


•iicr’VX 
m-<*  / 


.  Ar.e^  N 
nm  '<  ’ 


2  2 

applicable  In  the  case  v  »  o)  . 

These  graphs  can  be  used  to  determine,  from  the  phase  character¬ 
istic  of  the  atmospherics,  the  values  (N/v)gj.j.  vs.  frequency  (see 
below  and  [11]). 

Using  the  found  values  of  kS^  we  calculated  the  curves  of  the 
module  and  the  argument  of  Integrand  (l),  which  are  the  au^jlltude 
and  the  phase  characteristic  of  atmospherics  for  various  distances, 
viz.,  we  constructed  the  curves 


1  Q*  K  0)  (u,.  r)  I  =  B  («.,  T)  Ao  (u).  0)  =.  A  (<«.  r). 
•fg  (Qo  (“.  0)  u  (“>.  /■)!  =  ^0  (“*.  0)  —  '■)  ”  (<?  (“•  '■)• 
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Fig.  2.  Same  as  Fig.  1  for  frequencies 
of  15,  20,  and  30  kc  =  70,  60,  50  km), 
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Fig.  3.  Curves  of 

N/v  for  frequencies  of 
amd  10  kc  »  70  k 


Pig.  4.  Same  as  Pig.  3  for  frequencies  of  15,  20, 
and  30  kc  =  70,  60,  50  km). 


Pig.  5.  Curvea  of  |p(Cq)|  vs.  N/v  for 
frequencies  of  5»  7*  and  10  kc  = 


SB? 


Fig.  6.  Same  as  Fig,  5  for  frequencies 
of  15,  20,  and  30  kc  (h  =  70,  60, 

50  km). 


Fig,  7,  Curves  of  arg  P(C^j) 

vs.  N/v  for  frequencies  of  5, 
7,  and  10  kc  =  70  km). 
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Prom  the  flares  It  Is  evident  that  the  amplitude  curves  have 
two  maxima;  one  near  50-200  cps  euid  a  second  near  5-12  kc.  The  first 
of  these  Is  displaced^  with  distance,  toward  lower  frequencies}  the 
second  Is  displaced  toward  higher  frequencies.  The  first  maximum 
exceeds  the  second  by  a  factor  of  2-5,*  Near  2-5  kc  the  amplitude 
curves  have  a  minimum,  the  depth  of  which  Increases  greatly  with  dis¬ 
tance,  The  nature  of  the  pliase  curves  Is  also  quite  evident  from 
the  figures. 

The  number  of  waves  n  used  In  the  calculations  for  various  fre¬ 
quencies  and  distances  was  determined  such  that  their  modules  were 


*  Here  we  should  note  that,  as  follows  from  [1],  the  Ionosphere 
model  selected  for  the  calculations  Is  hardly  suitable  for  f  ^  3-4  kc, 
since  at  these  frequencies  It  Is  apparently  Impossible  to  tsdce  Into 
account  the  Infiltration  of  waves  through  the  Ionosphere.  Further, 
for  f  <  200-500  CDS  the  soxirce  spectrum  has  not  yet  been  Investigated, 
Therefore,  [g(<o,r;Q(co,0)|  and  arg[g((D,r)Q(aj,0) ]  calculated  here  require 
further  refinement  for  f  <  3-4  kc.  From  various  data  we  can  assxane 
that  the  first  maxima  In  the  frequency  range  50-200  cps  decrease 
several-fold  after  more  precise  calculations. 
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no  less  than  0,5^  of  the  modxae  of  the  most  Intense  vave  of  the  null 
node.  Although  In  certain  cases  It  was  equal  to  6-7  (for  r  =  500  kn), 
nevertheless  the  curves  in  Pigs.  9-12  are  reasonably  smooth,  since 
the  Increase  In  the  field  with  decreasing  frequency  smooths  the  com¬ 
plex  Interference  nature  of  the  change  In  field. 
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Pig.  11.  Same  as  Pig.  9,  for  r  =  2000  \m 
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Pig.  12.  Same  as  Pig,  9,  for  r  =  3000  ion. 

To  calculate  Integral  (l),  whose  Integrand  Is  given  using  Plg« 
9-12,  we  used  numerical  integration.  The  calculations  were  made  ir 
two  ways. 

In  one  of  these  we  used  the  method  of  harmonic  analysis  of  the 


Integrand;  this  method  Is  as  follows. 

Let  us  reirrlte  (7)  In  the  form 

£  (t,r)  =  jp  B  (u>.  r)  At  (u),  0)  cos  9  cos  Wdw  —  ( I5  ) 

"I 

—  i j  B  (u),  r)  Ag  (u),  0)  sin  <p  sin  widtu. 

“I 

If  we  now  determine  the  coefficients  aj^,  hj^,  Cj^,  and  dj^  of  the  Fourier 
series  of  the  expressions 


fl  (w,  r)i4g(u),0)cos(p  =  ^  I-  y  fl/icos  —  w  + 

^  Im 

.Vi.  ■ 

Sin  — 

h  “‘0 

B(c,r)  Ag(<o,0)  Sin  f  -  ^  i  ]£c*cos— 01+  yd*sin— w, 

k  V  “• 


(14) 


entering  Into  Integrals  (13),  we  can  obtain  the  following  formula  for 
calcvilatlng  E(tjr): 


c/,  V  •  I  .  V  —  COS‘“|.0  — a*sln“.< 

^  '■)  ■  *  L/  ^  ^ 

1  V  W„  / 

COS  «*>  /)  +  dk  tintaJ 

i 


(15) 


+ 


1  A  [  V'*!**  fc*  (1  —  cos«> /)  +  d*  1 

+  . - 1. 


which,  when  t  2TrVa)^,  assumes  the  form 


while  when  t  =  0 


£(/  =  0./-)=^<i,. 


(16) 


The  coefficients  for  expansion  Into  Povirler  series  can  be  deter¬ 
mined  by  haxTOonlc  analysis  of  curves  (l4). 

The  second  method  consisted  In  the  approximate  Integration  using 
the  Slii5)Son  fomula  by  breaking  down  the  Integration  Interval  Into 
the  required  number  of  parts.  For  this  the  Integrand  was  tabulated 
for  the  Interval  f  =  5-100,000  cps  for  various  values  of  t.  When 
using  the  Simpson  method  we  can  calculate  E(t,r)  for  any  values  of  t 


-16- 


(as  large  as  we  wish),  which  cannot  be  done,  for  all  practical  pur¬ 
poses,  using  the  first  method,  since  for  higher  harmonics  the  error 
In  the  harmonic  analysis  greatly  Increeuses. 

The  calculations  showed  that  In  the  range  f  =  1000-50,000  cps, 
which  corresponds  to  a  signal  duration  of  1  msec,  both  methods  give 
corresponding  results.  However,  to  calculate  the  so-called  slow 
"tails”  of  the  atmospherics  caused  by  the  propagation  of  waves  of 
superlow  frequency  f  «  500-1000  cps,  only  the  second  method  can  be 
used. 


2.  Calculation  Results  i  Comparison  With  Experimental  Data 

Using  the  method  described  In  the  previous  section  we  calculated 
the  signals  for  r  =  500,  1000,  2000,  and  3000  km  shown  In  Pigs.  13-16. 
The  flgiu?es  show  only  the  so-called  "high-frequency"  part  of  the 
atmospherics  (examined  In  this  work)  of  0.5-1  msec  In  duration  having 
the  fonn  of  damped  oscillations  with  Increasing  period.  The  number 
of  half-periods  of  this  part  of  the  signal  changes,  with  Increasing 
dlstauice,  from  appz*oxljnately  5  to  8}  the  duration  of  the  half-period 
Is  within  the  limits  of  30-200  msec.  Behind  the  high-frequency  part 
of  the  signal  there  often  follows  the  so-called  "low-frequency"  part, 
the  dTiratlon  of  which  often  reaches  10  msec  and  more  (see  Interval  b 
In  Fig.  17).  This  part  of  the  atmospherics  (the  so-called  "tall")  Is 
also  of  an  oscillatory  nature  but  with  a  much  lau?ger  half-period  of 
oscillations,  the  number  of  which  chamges  from  1  to  3  [12]  amd,  as 
cam  easily  be  \mderstood.  Is  formed  mainly  by  waves  of  frequency 
f  «  500-1000  cps. 

To  compare  the  calculated  signals  with  the  observed  atmospherics 
we  examined  504  oscillograms  recorded  from  September  1954  through 
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[illegible] 


Fig,  15.  Comparison  of  Pig.  14.  Same  as  Pig,  15, 

atmospherics  (a)  with  the  for  r  =  1000  km, 

retlcally  calculated  sig¬ 
nal  (b)  for  r  =  500  km. 


August  1956,  l.e,,  almost  two  years.  Atmospherics  having  Irregular 
form  were  not  Included  In  the  calculations;  we  examined  only  smooth 
atmospherics  of  type  A  [2],  Thus,  of  the  ^04  atmospherics  we  selected 
only  444  smooth  atmospherics.  Taking  Into  account  the  accuracy  In 
determining  the  distance  to  the  source  of  the  discharge,  the  examined 
atmospherics  were  grouped  in  the  following  distance  groups:  400-600, 
800-1200,  1800-2200,  and  2700-5500  km;  atmospherics  between  these 
selected  Intervals  were  not  examined. 


[Illegible] 


Pig.  15«  Same  as  Pig,  15, 
for  r  *  2000  km. 


Fig.  16.  Same  as  Fig.  15, 
for  r  =  5000  km. 


Plgixres  15-16  give,  for  comparison,  oscillograms  of  atmospherics 
recorded  at  approximately  these  distances.  Prom  the  figures  It  Is 
evident  that  there  Is  reasonably  good  external  slmll6u?ity  between  the 
theoretically  calculated  signals  and  the  atmospherics)  this  Is  expres¬ 
sed  In  the  closeness  of  the  number  of  transitions  n  through  the  zero 
point,  the  duration  of  the  high-frequency  part  of  the  atmospherics 
T,  and  the  attenuation  factor  a  (Pig.  18). 

A  con^jarlson  of  the  signals  according  to  these  external  character¬ 
istics  was  made  on  the  basis  of  all  experimental  data.  For  this  we 
constructed  the  distribution  curves  for  n  and  t  (Pigs.  19-22).  A 
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corresponding  analysis  for  the  value  of  a  was  made  only  for  the  most 
characteristic  days,  l.e.,  those  when  atmospherics  similar  to  the 
theoretical  signal  were  observed,  for  the  most  part. 


Fig.  17.  Photo-oscillogram 
of  atmospherics  with  the 
low-frequency  part. 


Fig.  18.  Form  of  the  atmos 
pherlcs  usually  observed  at 
r  >  2000-5000  km. 


Fig.  19.  Distribution  curves  of  the 
number  of  transitions  n  through  zero 
and  the  duration  t  of  ^e  atmospherics 


for  r  =  500  km. 


Fig.  20.  Same  as  Fig.  19, 
for  r  =  1000  km. 
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Fig.  21.  Same  as  Fig.  I9,  Fig.  22.  Same  as  Fig.  19» 

for  r  ==  2000  km.  for  r  =  30OO  km. 


The  distribution  of  the  examined  atmospherics  by  distances,  and 
also  their  number  In  per  cents,  which  was  similar  to  the  theoretical 
according  to  the  Indicated  characteristics  Is  given  In  Table  5. 

TABLE  3 


Oiatanoa 
Intorrala,  lai 

ggigglj 

Niabar  of  alallar 
anaa  froa  tha  nabbe 

88 

jjHMjll 

37 

74 

167 

21 

57 

1800-2200 

201 

200 

1 

48 

2700—3500 

48 

47 

1 

13' 

T0TM.I 

504 

444 

60 

*  Am  d*er«a««  In  tta*  nuriMr  of  nlallar  atwaphirlet  mX  r  "  2000-3000  la  la 
axplalmd.  In  partlouUr,  by  tka  fact  that  hart  «a  eonaldartd  atMapharlea  of  tha 
typa  Bhaan  In  Tie*  18  to  ba  dUaiallar  to  tha  thaoratloal  aiffialf,  althou^  tha 
eoiraspoiidlJig  aaluos  of  ft  end  t  ware  oloaa  (aaa  balow,  iftblo  8). 
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The  theoretically  calculated  signals  (Figs.  13-16)  have  the 
values  of  n,  o,  and  x  (Fig.  18)  given  In  Table  4.  From  Table  4  we 
see  that  n  and  t.  In  general,  are  little  sensitive  to  a  change  In  r. 


TABLE  4 


f.  KM 

,(  +  ) 

,(  +  ) 

500 

5 

7 

4,6 

10,8 

1000 

7 

6—7 

— 

0,8 

8,0 

2000 

7 

— 

~7 

— 

— 

7.0 

6,5 

3000 

8 

1 

4-5 

7.6 

— 

9,1 

,  13,2 

A  comparison  of  the  values  of  n  and  x  obtained  from  the  maxima 
of  the  distribution  curves  In  Figs.  19-22  with  the  data  In  Table  4 
gives  generally  satisfactory  agreement.  However,  beginning  at  dls- 
temces  r  ~  2000-3000  km  we  often  note  atmospherics  with  branches  which 
rise  and  fall  with  time  (Fig.  l8),  while  theoretically  such  a  fom  Is 
obtained  only  for  r  ~  3000  km,  and  rarely,  at  that.  However,  n  and 
T  are  close  to  their  theoretical  values. 

Naturally,  a  more  profound  analysis  of  the  similarity  of  the 
calculated  and  observed  signals  should  be  based  on  an  examination  of 
their  phase  and  an9>lltude  characteristics.  At  the  same  time,  since 
an  analysis  of  atmospherics  makes  It  possible  to  determine  only  the 
behavior  of  the  relative  amplitude  -5/"*'  which  changes  little 

with  distance.  It  Is  possible  to  obtain  new  data  only  from  curves  of 
the  additional  phase  ^((o,r),  which  changes  noticeably  with  distance. 
Since  from  the  experiments  It  Is  Inqjosslble  to  determine  the  absolute 
phase  values,  the  experimental  phase  curves  must  be  fitted  to  the 
theoretical  curves  In  the  region  of  maximum  frequencies  of  the  exam¬ 
ined  range,  where  changes  little  with  frequency  [1]. 
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As  Illustration,  Pig,  23  compares  the  behavior  of  the  phase 
characteristics  and  Pig.  24  the  behavior  of  the  amplitude  character¬ 
istics  of  two  externally  similar  atmospherics  at  distances  of  approxi¬ 
mately  1000  and  3000  km.  Prom  Pig,  23  It  Is  evident  that  the  phase 
curves  of  the  experimental  atmospherics,  fitted  to  the  theoretical 
curves  at  a  frequency  of  10-15  kc,  are  generally  close  to  the  theo¬ 
retical  curves  right  up  to  3  kc.  We  must  mention  that  In  the  majority 
of  cases,  at  frequency  f  <  3  kc  there  Is  divergence  of  the  phase  and 
anq>lltude  curves.  All  this  attests  to  the  fact  that  the  effective 
Ionospheric  parameters  selected  at  these  frequencies  evidently  do  not 
correspond  to  the  parameters  of  the  real  Ionosphere. 


Pig,  23,  Conqjarison  of  the  theoreti¬ 
cal  (solid  lines)  and  experimental 
(dashed  lines)  amplitude  characteris¬ 
tics. 
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A  comparison  of  the  phase 
characteristics  of  the  theoreti¬ 
cal  and  experimental  atmospherics 
Is  of  great  Interest  In  that,  as 
Al'pert  has  shown  [1],  It  makes 
Biu.r)A,(w.o}  It  possible  to  determine  the  effec 

tlve  parameters  of  the  Ionosphere 
and  the  rate  of  propagation  of 
low-frequency  waves.  Since  fioifr) 
Is  a  function  of  the  Ionospheric 
conductivity,  proportional  to  N/v, 

Pig.  24.  Con^jarlson  of  the  theo-  where  N  Is  the  electron  concentra- 
retlcal  (solid  lines]  and  e3q>erl- 

mental  (dashed  lines)  an^lltude  tlon  and  v  Is  the  number  of  col- 
character  Is  tics  . 

llslons,  with  certain  assumptions 
we  can,  from  the  curves  V'(a>»r),  determine  the  values  of  Sq^^  (Pigs.  3- 
4)  deteimiine  N/v  from  data  on  the  analysis  of  one  or  a  group  of  atmos¬ 
pherics  recorded  \mder  any  specific  conditions.  With  a  high-speed 
harmonic  analyzer  we  can  obtain  these  data  quite  rapidly,  l;e.,  use 
the  corresponding  method  for  practical  purposes.  This  question  has 
been  examined  In  more  detail  by  Al»pert  [l]|  here,  however.  It  Is 
Interesting  to  make  a  more  detailed  analysis  of  the  Individual  atmos¬ 
pherics  received  on  various  days. 

We  noted  cases  when  throughout  the  day,  practically  from  the 
same  direction  and  distance,  we  received  atmospherics  of  Identical 
form,  often  similar  to  the  theoretical  signal.  This  made  possible 
the  conclusion  that  throughout  the  day  the  effective  parameters  of 
the  Ionosphere  along  a  given  atmospherics  propagation  path  were  rea¬ 
sonably  constant  and  close  to  those  obtained  In  the  calculations. 
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For  example,  on  July  18,  1955#  we  recorded  10  atmospherics,  from  dis¬ 
tances  of  400-600  km,  from  0950  hours  to  1750  hours.  Two  of  these 
were  of  Irregular  form  and  were  not  examlnedi  the  remaining  eight 
were  all  of  Identical  form.  It  Is  characteristic  that  all  atmospher¬ 
ics  except  the  last  came  from  the  same  direction  and,  consequently, 
had  the  same  thunderstorm  center  as  the  soxirce.  In  Table  5  we  give 
the  parameters  of  the  atmospherics  recorded  during  days  when  all  smooth 
atmospherics  were  generally  of  approximately  the  same  form  as  the 
*■  theoretical  signals. 


TABLE  5 


OiTM- 

Ota- 

H 

-(  +  ) 

•mw 

Itet« 

MSfetlM 

tlM 

taiaa. 

Imaln 

kail 

>«•  aao'* 

io*aae 

/■  ~  500  (Ml 

18.VI.18SS 

9  49 

a 

500 

7 

0 

4,3 

3,2 

12  46 

s 

600 

7 

7 

5,6 

5,0 

12  SO 

s 

600 

7 

7 

7.3 

4,8 

12  SI 

s 

500 

6 

7 

5,7 

9,7 

12  S4 

sa 

500 

7 

7 

4,6 

5,6 

17  46 

s 

600 

6 

6 

6,3 

3,6 

17  46 

s 

500 

5 

5 

8.4 

7.2 

17  S3 

s 

600 

S 

5 

8,8 

6.4 

10.VI1I.1956 

12  49 

s 

500 

5 

15 

1,3 

1.5 

16  4S 

NU 

400 

8 

12,5 

6.1 

4,2 

16  49 

MNH 

600 

8 

10 

4,5 

5.2 

30.VIII.19S6 

9  SI 

s 

450 

6 

9,5 

3,6 

4.1 

'  9  54 

1  5  i 

450 

7 

10,5 

5.0 

7,9 

J 

r  ~  1000  KJM 

24.Vl.19SS 

17  S2 

SSE 

10.50 

6 

4.4 

5,0 

8,5 

17  56 

SMbM 

975 

9 

7 

6,5 

5,4 

22.VI.19SS 

12  S3 

SSN 

1100 

6 

5 

8,6 

5,6 

17  45 

SSH 

900 

6 

4,5 

6.0 

11,8 

17  51 

i  s 

900 

6 

4,5 

5,8 

5.2 

17  54 

s 

975 

5 

4 

13,7 

5.6 

8.VI.19S6 

9  46 

SUM 

1100 

7 

8,5 

9  SO 

SUM 

1100 

9 

11 

— 

» 

12  45 

S 

7.50 

8 

12.5 

3,6 

4,4 

26.VI.19S6 

9  46 

ME 

750 

7 

10 

9  47 

NE 

850 

6 

5,5 

•  4,9 

7.2 

9  48 

ME 

7.50 

6 

6,5 

6,4 

7,5 

9  50 

S 

1100 

5 

3.9 

6,0 

2,9 

9  52 

S 

1000 

7 

7.5 

7,2 

7.9 

4.VII.19S0 

12  49 

H 

1100 

6 

5,5 

4.6 

6.3 

16  48 

U 

1100 

7 

8 

3.6 

5.5 

16  48 

ssw 

8.50 

6 

7 

4,0 

2.8 

16  S3 

H 

1100 

8 

11,5' 

2.4 

4.5 

16  54 

w 

1050 

8 

10,5 

5.2 

4.7 

25- 


TABLE  5  (continued) 


1 

Date  1 

DUree- 

tlon 

n 

tf. 

10^  SM 

,(+) 

io*.eao-> 

io«  eee-* 

r 

~  2000  KM 

27.IX.1954 

9  50 

SSE 

1800 

I  ^ 

1  2,2 

1,3 

— 

9  54 

SSE 

2000 

3 

10,6 

8.3 

10.11.1954 

12  53 

ssu 

2200 

6 

3.5 

6.3 

3,8 

12  55 

SSH 

1800 

7 

5 

5,6 

571 

17  53 

SU 

1800 

6 

4.5 

5.4 

7.0 

12.IV.1955 

9  47 

SW 

2200 

8 

6 

4.1 

4,8 

5  50 

ssw 

2000 

8 

9 

3.2 

3,0 

9.V.105S 

12  49 

c 

1950 

6 

3,5 

5,8 

7.1 

17  45 

ssw 

1870 

6 

4 

7.7 

10,8 

22.X.1955 

9  46 

s 

2000 

7 

5,5 

6.1 

7,7 

9  49 

s 

2100 

5 

3.5 

16.7 

9,2 

< 

12  46 

s 

2100 

7 

— 

— 

— 

16  48 

s 

1900 

6 

4 

6.0 

6,1 

23.111.1956 

9  49 

SSK 

itoo 

7 

3.7 

3,4 

10,0 

9  49 

SSI 

1875 

5 

2,5 

20,7 

9.4 

9  55 

s 

2025 

8 

6 

6,5 

8.4 

12  52 

SSI 

1950 

6 

3.3 

6,8 

3,3 

r  ~  3030  KM 

18.IX.1954 

12  50 

t 

3500 

8 

2.0 

8,3 

6.2 

12  57 

KlfeC 

3300 

8 

4.3 

5,5 

11,9 

Comparison  of  the  data  In  Tables  3  and  4  shows  quite  good  agree¬ 
ment  between  n,  t,  and 

When  Identical  atmospherics,  similar  to  the  theoretical  signal, 
were  recorded  at  different  times  and  from  different  distances  and 
directions,  the  effective  Ionospheric  parameters  were  evidently  Iden¬ 
tical  for  a  long  period  of  time,  not  only  along  one  path  but  In  a 
rather  broad  region.  For  example,  such  was  the  case  on  Jvine  23,  1936, 
at  distances  of  300  and  1000  km.  On  this  day,  no  other  atmospherics 
were  observed  (irregular  form,  nonslmllar)  (Table  6). 

TABLE  6 


Date 

Meaaure- 

Hnt  tlM 
Iriatn 

Oiree- 

tlen 

Die- 

tanee, 

kM 

n 

^1. 

I0-*  Me 

,(  +  ) 

tar 
io>  aac> 

.<-> 

ter 

lo-eer* 

25.  VI.  1956 

12  S3 

I 

500 

6 

7 

5,2 

11,3 

16  51 

SM 

lo:x) 

7 

8 

8.2 

3,11 

16  52 

E 

375 

5 

5 

5,8 

8,4 

16  54 

SE 

400 

5 

5 

8,9 

8,4 

16  56 

s 

850 

G 

6 

6,3 

4,7 
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Another  such  case  Is  Illustrated  In  Table  7.  On  this  day  one 


atmospheric  of  Irregular  form  was  also  recorded. 


TABLE  7 


Date 

r*'lWV^ 

Dla- 

It 

't. 

,(  +  ) 

i(4ao^* 

icVK-* 

22.  VI.  1955 

12  47 

sw 

1870 

8 

5;0 

i 

14,7 

5,5 

12  S3 

ssw 

1100 

6 

4.0 

3.3 

5,6 

17  45 

ssw 

900 

;6 

4.5 

;o,o 

11,8 

17  49 

SH 

1870 

,7 

5 

!5.7 

4.8 

17  51 

s 

900 

,6 

4.5 

5.8 

5.2 

17 

sw 

1800 

,6 

5 

5.0 

7.3 

17  54 

t 

s 

.  975 

'5 

1 

4, 

1(3.7 

5,6 

On  certain  days  all  the  recorded  atmospherics,  or  peui:  of  them, 
differ  externally  from  the  calculated  ones,  and  dissimilar  atmospher¬ 
ics  were  often  received  from  the  same  direction  and  distance,  and  at 
practically  the  same  time.  For  exan^le,  on  Jime  ^0,  1956,  six  atmos¬ 
pherics  from  a  distance  of  approximately  1000  km  were  recorded.  Five 
of  these  were  similar  to  the  theoretical  signal  and  one  differed  from 
It,  although  they  were  all  recorded  at  approximately  the  same  time 
and  from  the  same  sotirce.  Another  such  case  occurred  on  July  21,  1956. 

At  present  the  reasons  for  such  cases  are  not  clear. 

Table  8  gives  the  distribution,  by  days,  of  similar  and  dissimilar 
atmospherics.  The  second  column  gives  the  total  number  of  days  when 
atmospherics  were  observed  from  the  given  distance.  The  third  column 
shows  the  number  of  days  when  only  smooth  atmospherics  were  observed 
(since  there  were  days  when  only  Irregvilar  atmospherics  were  observed, 
and  these  were  not  Included  In  the  examination).  The  fourth,  fifth, 
and  sixth  columns  give  the  nvanber  of  days.  In  per  cents,  when  we 
observed  only  similar  atmospherics,*  similar  and  dissimilar  atmospherics. 

The  number  of  similar  atmospherics  In  Table  8  also  Include 
atmospherics  of  the  type  shown  In  Fig.  l8.  If  their  n  and  t  values  are 
close  to  the  values  of  n  eind  t  for  the  theoretically  calculated  atmos¬ 
pherics. 
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l.e.j  mixed,  and  only  dissimilar  atmospherics.  From  Table  8  It  Is 
evident  that  the  niimber  of  days  when  similar  atmospherics  were  observed 
decreases  with  Increasing  distance.  As  we  have  already  mentioned, 
the  form  of  the  atmospherics  noted  at  distances  of  2000-3000  km  are 
obtained  theoretically  only  for  a  distance  of  3000  km.  This  can  pos¬ 
sibly  be  explained  by  the  fact  that  In  the  ceCLculatlons  not  all  the 
features  of  propagation  which  begin  to  appear  at  great  distances  were 
taken  Into  consideration. 


TABLE  8 


Dlatano* 

interval* 

\m 


4CW— COO 
800—1200 
1800—2200 
ZiOO-liriOO 


Mt«l  BU^ 
^r  of  days 
Uten  atsMs- 
^harles  wti 
obsarvad 


iiwl 


MUBbar  of 
M[^s  ahan 
smooth  atiaoa^ 
iherles  wora 

obsarvad 


Nisabar  of 
days  when 
only  slal* 
Isr  atmos- 
pharlos  aa 


Nupaber  of 
dsys  i&mn 
■ixad  atne» 
pharles  aer|, 
obtyvad] 


Muidiar  of 

isitsn 

liar  stMs> 
phonos 

sbsarvad.  SI 


Hunbar  of 
atnosphei^ 
las  with 
two  similar 
oharmotprlstlo* 
(jn  aadt),  S 


38 

(Kl 

H(> 

31 


26 

57 

80 

30 


61,5 

45 

31 

7 


11,5 

2:1 

31 

10 


27 

32 

37 

83 


80 

01 

60 

57 


Conclusion 


The  results  of  our  work  show  that  the  signals  calculated  theo¬ 
retically  are  In  good  agreement  according  to  all  external  character¬ 
istics.  In  40-60J^  of  the  cases,  the  similarity  was  not  only  quali¬ 
tative  but  also  quantitative,  which  Indicates  that  In  these  cases  the 
effective  Ionospheric  parameters  selected  for  the  calculations  were 
apparently  realized.  This,  in  turn,  confirms  the  earlier  assunqjtion 
[l]  that  as  a  result  of  the  analysis  of  the  phase  characteristics  of 
atmospherics  we  can  determine  the  effective  conductivity  of  the  Ionos¬ 
phere  during  the  day  for  specific  conditions.*  There  Is  a  certain 

■  '*  oJhe  present  work  Is  limited  to  ein  examination  of  atmospherics 
observed  d\irlng  the  day,  since  after  sunset  the  field  Is  formed.  In 
the  majority  of  cases,  by  several  modes,  and  it  becomes  difficult  to 
compare  theoretical  data  with  atmospherics  analysis  results. 
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discrepancy  between  the  form  of  the  atmospherics  and  the  theoretically 
calculated  signals  with  Increasing  distance,  which  warrants  further 
analysis. 

In  conclusion  I  should  like  to  express  my  thanks  to  Ya.  L. 

Al’pert  for  suggesting  this  question,  for  his  leadership,  and  for 
his  constant  help  with  my  work,  to  S.  V,  Borodina  for  her  advice 
during  my  work,  and  to  V.  I.  Krayushklna  for  her  help  with  the  cal¬ 
culations  and  setting  up  the  tables  and  graphs. 

PTCWERENCES 


1.  Ya.  L.  M*pert.  Radlotekhnlka  1  Elektronlka,  1,  No. 

1956;  S.  V.  Sorodlna,  Zhumal  Eksperimental •  noy  1  Teoretlcheskoy 
Plzlkl,'  1957:  Radlotekhnlka  1  Elektronlka,  4,  195,  1959. 

2.  S.  V.  Borodina.  Trudy  Nauchno-Issledovatel'skly  Instltut 
Zesmogo  Magnet Izma,  lonosfery,  1  Rasprostranenlya  Radlovoln,  No.  i;5, 

1956. 

3.  Ya.  L.  iU.*pert.  The  Propagation  of  Low-Prequency  Electro¬ 
magnetic  Naves  Above  the  Earth* s  Surface.  Moscow,  Izdatel’stvo 
Akademll  Nauk  SSSR,  1955. 

4.  C.  Bruce,  and  R.  Oolde.  J.I.E.E.,  88,  48?,  1941. 

5.  J.  Barlow,  Q.  Prey,  and  J,  Newman.  J.  Pranklln  Inst.,  258, 

187,  1954. 

6.  K.  Budden.  Phil.  Mag.,  42,  I-I9,  1951. 

7.  J.  Walt.  Natl.  Bur.  Stan.  Rept.  5515,  September  3,  1957. 

8.  Ya.  L.  Al»pert,  and  D.  S.  Pllgel*.  Radlotekhnlka  1  Elek¬ 
tronlka,  4,  202,  I959. 

9.  Ya.  L.  Al*pert.  Uspekhl  Plzichesklkh  Nauk,  60,  369,  I956. 

10.  Table  of  the  Bessel  Punctlons  Jq(z)  and  J2^(z)  for  Con5)lex 

Arguments.  Columbia  University  Press,  1943;  Table  of  the  Bessel 
Punctlons  Y(z)  and  Y(z)*  for  Complex  arguments.  Columbia  University 
Press,  1950. 


*  As  In  original  [Translator's  note]. 


PTD-TT-62-1193-1+2+4 


-29- 


11.  L.  M.  Brekhovsklkh.  Doctoral  Dissertation.  Plzlcheskly 
Instltut  Alcademll  Nauk  SSSR«  1946;  Izvestlya  Akademll  Nauk  SSSR, 
Serlya  Plzlcheskaya,  8,  505;  15,  515,  594,  1949. 

12.  L.  Lleterman.  J.  Appl.  Phys.,  27,  1475,  1477,  1956. 


PTD-TT-62-1195-1+2+4 


-50- 


DISTRIBUTION  LIST 


DSPARilfflNT  OP  L3PSNSE  Nr.  Copies 


iiSADQUiVRTERS  USAP 
apci:j-3b2 

(ARB) 


OTlifiJR  AGENCIES 


C  lA 

AID 

0T3 

'  \7  P 

P'>.'3 

NASA 

RAND 


1 

6 

2 

2 

2 

1 

1 

1 


LiiUOR  AIR  COlffiLAIIDS  Nr.  Copie 
APSC 

SCPTR  I 

ASTIA  25 

TD-BIa  5 

TD-Blb  3 


ADBC  (AEY) 
SSD  (SSP) 
APGC  (PGP) 
BSD  (ESY) 
RADC  (RAY) 
APSV/C  (SWP) 
APMTC  (mTW) 


62-1193/1+2+4 


31 


